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ABSTRACT 

The results of spectral modeling of the data for a series of RXTE observations and 
four ASCA observations of GRO J 1655-40 are presented. The thermal Comptonization 
model is used instead of the power-law model for the hard component of the 
two-component continuum spectra. The previously reported dramatic variations of the 
apparent inner disk radius of GRO J1655-40 during its outburst may be due to the 
inverse Compton scattering in the hot corona. A procedure is developed for making 
the radiative transfer correction to the fitting parameters from RXTE data and a more 
stable inner disk radius is obtained. 

A practical process of determining the color correction (hardening) factor from 
observational data is proposed and applied to the four ASCA observations of 
GRO J 1655-40. We found that the color correction factor may vary significantly 
between different observations and the finally corrected physical inner disk radius 
remains reasonably stable over a large range of luminosity and spectral states. 

Subject headings: accretion, accretion disks — black hole physics — stars: individual 
(GRO J 1655-40) — X-rays: stars 


1. Introduction 

The X-Ray Nova GRO J1655-40 (Nova Scorpii 1994) was discovered by the Burst and 
Transient Source Experiment (BATSE) onboard the Compton Gamma Ray Observatory on July 
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27, 1994 (Zhang et al, 1994). Superluminal jets were soon observed at radio frequencies from the 
system (Tingay et al., 1995; Hjellming and Rupen, 1995) in close correlation with its hard X-ray 
activities (Harmon et al., 1995), thus suggesting its similar nature to the black hole candidate 
and the first microquasar GRS 1915+105 (Mirabel and Rodriguze, 1998). Optical observations 
allowed precise determination of the parameters of the binary system: mass of the compact object 
of 7.02 ± 0.22 M q , inclination of i = 69.50 ± 0.08° degrees and a mass ratio of Q = 2.99 ± 0.08 
(Orosz & Bailyn 1997). Similar, yet more conservative estimate of the black hole mass obtained 
by Shahbaz et a/.(1999) is 6.0 ± 0.2 M 0 . These results thus provide so far the best evidence for 
the black hole nature of the compact object in GRO J1655-40. 

The proximity of GRO J 1655-40 and its well determined system parameters make it one of 
the best laboratories for studying several important astrophysical subjects: accretion disk physics, 
jet/outflow production and transportation, general relativity in the strong field limits, etc. Some 
of these results have already been reported. Hameury et o/.(1997) suggested that the first observed 
optical brightening prior to an X-ray outburst in GRO J1655-40 (Orosz et al., 1997b) provides 
evidence for the advection dominated accretion flow (ADAF) in the quiescent state of black hole 
X-ray binaries (Narayan and Yi, 1995; Narayan, 1995). Zhang, Cui and Chen (1997b) made the 
first determination of the black hole angular momentum of GRO J1655-40, making use of the high 
quality energy spectrum of the source observed with ASCA (Zhang et al., 1997a). The estimated 
near maximal black hole spin rate in GRO J1655-40 and the other microquasar GRS 1915+105 
suggests that the highly relativistic jets may be related to the angular momentum of the black 
holes (Blandford and Znajek, 1977). Cui, Zhang and Chen (1998) have also proposed that the 
~300 Hz QPO observed in GRO J1655-40 (the highest frequency from any black hole system) is 
the frame dragging effects around a rapidly spinning black hole, predicted in Einstein’s general 
relativity. If confirmed, this would be one of the first cases of observational verification of general 
relativity effects near the strong field limit. More recently, Zhang et al. (2000) discovered the 
three-layered accretion disk structure in the two microquasars GRO J1655-40 and GRS1915+105, 
making use of the ASCA and RXTE data. The similarity between the accretion disk structure and 
the solar atmosphere suggests that magnetic fields may play very important roles in the accretion 
process, high energy electron acceleration, high energy radiation and outflow /jet production. 

However, one fundamental and unsolved question in the accretion disk physics is whether the 
inner disk boundary is stable, especially in the high luminosity state of black hole X-ray binaries. 
In the quiescent state the inner portion of the geometrically thin and optically thick disk is most 
likely truncated very far away from the black hole horizon, and replaced by a hot, optically thin 
and nearly spherical accretion flow. This is the well known ADAF solution. Alternative outflow 
may dominate inward accretion and the near-Keplerian disk may still extend very close to the 
black hole horizon (Blandford and Begelman 1999). Observationally the two scenarios are still 
indistinguishable, due to the limitations of the currently available X-ray instruments. In the high 
accretion state the disk is generally believed to be able to extend very close to the black hole 
horizon, because strong X-ray emissions are observed, indicating that the conversion efficiency 


from the gravitational potential energy into radiation energy is high, and thus the inner disk 
radius should not be very large. It is however not clear if the inner disk boundary may stay stably 
at the last stable orbit of the black hole, regardless of the luminosity and spectral states of the 
source. This issue is very important since if the inner disk radius varies significantly, one must be 
very careful in determining the radius of the last stable orbit, and thus the previously reported 
determination of black hole angular momentum may be questionable. 

Significant inner disk radius variations have been reported in GRS 1915+105 (Belloni et al., 
1997) and the variations are found to be consistent to the thermal and viscous time scales of 
the disk at various radii. However, recently Merloni et a/. (1999) suggested that these variations 
might be due to the variations of the disk structure, which resulted in the variations of the 
color hardening factor, and thus made the apparent and uncorrected inner disk radius highly 
variable. In GRO J 1655-40 and several other black hole X-ray binaries, the apparent inner disk 
radii inferred from X-ray spectral modeling have also been found to vary significantly. Here in 
this paper we report our re-analysis of some of these data, with the incorporation of a radiative 
transfer correction procedure, because some of the photons emitted from the accretion disk could 
not reach the observer directly; they are converted to hard X-ray photons by the hot electrons in 
the corona around the black hole. The details of the radiative transfer correction are presented 
in Appendix A. We also propose a procedure for the determination of the color hardening factor 
from data; the procedure is presented in Appendix B. Our results demonstrate that the inner disk 
radius in GRO J 1655-40 is consistent with being a constant; the apparent and uncorrected inner 
disk radius variations are due to both the radiative transfer and color hardening processes in the 
inner disk region. The angular momentum of the black hole in GRO J 1655-40 is found to be near 
maximal, consistent with early reports by some of us (Zhang et al . , 1997; Zhang, Cui and Chen 
1997; Cui, Zhang and Chen 1997). 


2. Data Selection and Analysis 
2.1. Data Selection 

Like typical BHC (black hole candidate) in X-ray binary system, GRO J1655-40 underwent a 
rising phase, a very-high/hard state, a high/soft state and a low/hard state in the 1997 outburst 
(Sobczak et al., 1999), as shown in Fig. 1. Some typical spectra in these states are shown in 
Fig. 2. We present 38 observations, as shown in Fig. 1 , covering primarily the very-high state (in 
which the source underwent a complex spectra variation) of the 1997 outburst of GRO J1655-40 
obtained with PCA (Proportional Counting Array) onboard RXTE; Sobczak et a/.(1999) reported 
significant inner disk radius variation during this period. The “standard2” mode data are used, 
which correspond to 129 channels from 0-100 keV. The PCA is made of five PCUs. Considering 
instrument response, only data from PCU0 and PCU1 in the energy range 2.5-20 keV are 
used (Sobczak et al., 1999). The response matrix is Keith Jahoda’s 1998 January version. The 


background spectra are obtained using the standard background model for bright sources. The 
spectra are extracted separately for each PCU and fitted simultaneously using XSPEC. A 1% 
systematic error was included in the fitting. 

Four observations of the bright state of GEO J1655-40 have been made with ASCA since its 
launch. Its excellent response in low energy range provides us with the possibility of studying the 
spectrum in this area. The GIS2 events were collected within a circular region with radii between 
5’ and 6’ centered at the source positions. Background was not substracted. The GIS2 data in 
the energy range of 0.7 - 10 keV are used in our analysis. For one of these observations we also 
have a simultaneous RXTE observation, which allows us to evaluate and compare the relative 
normalization of the RXTE/PCA and the ASCA/GIS instruments. This observation of ASCA 
lasted for two days. When doing the ASCA-XTE joint fitting, we only used the data collected 
simultaneously with RXTE. 


2.2. Fitting Models and Parameters 

Traditionally the continuum spectra of black hole X-ray binaries have been fitted with a 
two-component model, a soft component approximated by a disk-blackbody-like spectrum, and a 
hard component approximated by a power-law-like model with a high energy cutoff. For reasons 
discussed in Appendix A, we replace the power-law component with the thermal Comptonization 
model of Titarchuk (compTT in the XSPEC package, Titarchuk et al, 1995). In this model, the 
seed photons to the Comptonization process axe assumed to have a characteristic energy, which 
corresponds roughly to the lower end cutoff of the Comptonized component. The higher end cutoff 
corresponds to the characteristic energy of the hot electrons. The comparison of the two models is 
shown in Fig. 3. 

Thus the general model we employed for the continuum spectra is: wabs(diskbb + compTT). 
For the observations only with the RXTE instrument, the neutral hydrogen column density is 
fixed at 0.66, which is inferred from the ASCA and RXTE joint fitting of the simultaneous RXTE 
and ASCA observations, which occurred in the high/soft state of GRO J1655-40, in which the 
hard component can be ignored below 10 keV. In the rising phase and the high/soft state, an 
absorption edge is required for the fitting. Sometimes a line feature is also required . Because 
our main objective is to quantify the continuum spectra, the edge and line parameters are not 
presented in this paper. 

In the RXTE/ASCA joint fitting, for the two data sets we use the same model and parameters, 
except for the different scaling factors, i.e., the model becomes wabs(diskbb + compTT) constant, 
in which constant is fixed to unity for the RXTE data and allowed to float for the ASCA data. 
For three of the four ASCA observations, the soft component can also be fitted with the compTT 
model. The parameters of the two models for the same data will give us some understanding of 
the hardening process in the disk. 


The fitting parameters for the RXTE observations are listed in Table 5 and shown in Fig. 4. 
It can be seen that in the very-high/ variable state, the soft component normalization varies 
dramatically, and there is a significant anti-correlation between the hard photon flux and the soft 
component normalization. The fitting parameters for the ASCA observations are listed in Table 5 
and shown in Fig. 6. The ASCA-RXTE joint fitting results are shown in Table 2 and shown in 
Fig. 5. 


3. Results 

3.1. Interpretation of the multi-color disk black-body 

The inner radius of the disk around the black hole may be estimated if the parameters of 
the diskbb component are determined from the X-ray spectral fitting. In the soft state in which 
the hard component is negligible, we might consider that no photons emitted from the disk are 
scattered to the hard component, i.e., all disk photons escape from the system and reach the 
observer (albeit to the interstellar absorption and the limited detector area and response, which 
may be corrected after the observations). Thus the multicolor disk normalization parameter is 
related to the inner disk radius in the following equation (Makishima et ai, 1986): 



where /?;□ is the inner disk radius in kilometers, D is the distance to the source in kpc, 9 is the 
inclination angle of the disk. The geometric parameters of GRO J 1655-40 are well-determined 
from optical observations: 9 — 69.5° ± 0.1° (Orosz & Bailyn, 1997), D = 3.2 ± 0.2 kpc (Hjellming 
& Rupen, 1995). The inner disk radius can thus be determined from this relation. 

If the hard component is comparable or stronger than the soft component, the observed 
soft component cannot be considered as the original disk emission spectrum. This is because 
some of the disk photons are Compton scattered to form the hard component, and consequently 
the observed soft component is only a portion of the original disk spectrum. Therefore the 
normalization parameter inferred from the direct spectral fitting under-estimates the original disk 
emission significantly. Assuming that the disk photons are the main source of the seed photons of 
the Comptonization process responsible for the hard X-ray component, it is possible to make a 
radiative transfer correction in order to restore the original disk emission spectrum. The commonly 
observed hard and soft flux anti-correction provides a strong support to the assumption that the 
disk photons are the main source of seed photons. We have thus developed a procedure for the 
radiative transfer correction, as presented in Appendix A. Reliable hard component parameters 
are required for the radiative transfer correction, which can be obtained from the RXTE data. 
The ASCA data are usually not sufficient for determining the hard component reliably, because 
the ASCA instruments lack any useful response above 10 keV. 


However, an appropriate color correction should be also applied to get the physical inner disk 
radius, due to the saturated electron scattering of the original blackbody photons in the inner 
disk region, where electron scattering dominates the free-free absorption completely. Assuming 
that the electron scattering changes the temperature Teg of the soft component by a factor of /, 
i.e., the temperature becomes fT e g, the actual inner disk radius should be Reg = f 2 R\ B . The 
commonly used color correction factor from Shimura & Takahara (1995) is / = 1.7 ±0.2. However, 
the color factor might not be constant; it should change with the disk structure. Recently Merloni 
et al. (1999) suggested that the color correction factor might range from 1.5 to 3, thus indicating 
the estimated inner disk radius may be quite uncertain, unless the color hardening factor is 
known reliably. In Appendix B, we propose for a procedure to determine the color correction 
factor reliably from the same spectral data, provided that the instrument has sufficient low energy 
response and resolution, such as the ASCA instruments, but not the RXTE instruments. 

Thus the appropriate procedure for estimating the physical inner disk radius after the 
two-component spectral fitting should be: 1) do Comptonization correction if the hard component 
is not negligible, and 2) do color correction, if the color correction factor may be determined from 
the data reliably. 


3.2. Radiative Transfer Correction 

Fig. 4 and Table 5 shows the spectral fitting results for the 38 RXTE/PCA observations, as 
well as the corrected normalization parameter of the disk-blackbody component after the radiative 
transfer correction procedure is applied. It is clear that before the correction, the normalization 
parameter varies significantly and is anti-correlated with the hard X-ray photon flux. After the 
radiation transfer correction, the new normalization parameter, which is proportional to R 2 , does 
not vary very much, and thus implies that the inner accretion disk radius is probably stable. 
However, the inner disk radius cannot be estimated reliably from the PCA/RXTE data alone, 
because the precise value of the color correction factor cannot be obtained. 


3.3. Color Correction 

Though it is difficult to obtain hard component parameters reliably with the ASCA data, 
fortunately in all four ASCA observations the hard component is very weak, as determined from 
the simultaneous CGRO/BATSE and RXTE observations. Therefore the Compton radiative 
transfer correction is not important; only one of them needs minor radiative transfer correction. 
Fig. 6 shows the energy spectra from several ASCA observations. Table 5 shows the fitting 
parameters of the ASCA data, the value of color correction factor, the normalization parameters 
from these observations (see Appendix B for details). The data from the first observation cannot 
be fitted with the diskbb model, because the optical depth is small, and the energy difference 


between the incident photons (T 0 ) and the hot electrons (kT) is large. This suggests that the 
spectrum shape cannot be well described by a thermal spectral model, and the color correction 
factor is not meaningful in this case. It can be seen that the color correction factor is not a 
constant; the range of the color correction variation is consistent to the prediction of Merloni 
et al (1999). Because the exact value of the interstellar absorption is unknown, the estimated 
parameters suffer significant uncertainties, as shown in the table. 

Assuming that the inner disk radius is the last stable orbit of the black hole and following 
the procedure in Zhang, Cui and Chen (1996 Which is this?), we estimate that the black hole 
angular momentum of GRO 1655-40 is around between 0.8-1. 0 GAf/c 2 , consistent to an extremal 
Kerr black hole. 


3.4. Relative Calibration Issue of ASCA and RXTE 

We want to address briefly the relative calibration issue between RXTE/PCA and 
ASCA/GIS. Fig. 2 and Table 2 show the model-fitting results for the simultaneous RXTE 
and ASCA observations on July 26, 1997. Clearly there is a significant difference between the 
calibration of ASCA/GIS and RXTE/PCA for the soft component. For this observation the 
RXTE/PCA overestimated the low energy flux significantly, if we trust the absolute calibration 
of the low energy response (around 1 keV) of ASCA. One therefore must be very careful in 
interpreting the model fitting results using any RXTE/PCA low energy data (roughly below 5 
keV). 


4. Discussion 

The application of the thermal Comptonization spectral model compTT to the hard 
component and the subsequent radiative transfer correction procedure are only valid to a certain 
extend. If the source of the seed photons is not the soft component, the whole radiative transfer 
correction procedure is not meaningful at all. However, as shown in detail by Zhang et a/.(2000), 
the characteristic photon energy of the seed photons is very close to that of the disk photons, 
as determined with the ASCA data. This thus suggests that the dominant source of the seed 
photons is indeed the soft component, at least during these ASCA observations. However, if there 
is evidence for other significant source of seed photons, we must take this into account in the 
radiative transfer correction process. 

It has been shown from some high energy observations extending beyond several hundred keV 
that the high energy electrons responsible for the Comptonization process in the corona may have 
a non-thermal energy distribution. In this case the parameters derived for the compTT model 
are not meaningful physically. However, because the radiative transfer correction procedure only 
uses the total photons in the hard component, the final normalization parameter derived after the 
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correction process should still be valid. 

fundamental limitations of the radiative transfer correction procedure are related to the 
assumptions made for the properties of the corona, which is assumed to have a spherical shape 
with a density profile of l/R and a truncation outer radius constant during all observations. We 
have not studied quantitatively the uncertainties introduced by any deviations of the assumptions 
from the real corona. However, our results imply that the correction procedure probably has taken 
the main physical effects into account, because the corrected inner disk radius appears to be quite 
stable for most of the observations. One way for calibrating our model is to compare the corrected 
inner disk radius with that inferred from the soft state, in which no radiative transfer correction 
is needed. It can be seen from our data that by choosing an appropriate outer truncation radius 
of the corona, the corrected inner disk radius from all observations in the very-high/ variable state 
may be quite consistent with that inferred in the soft state (the inner disk radius in the soft state 
has always been found to be very stable). 

The procedure for determining the color correction factor is also subjected to a number of 
uncertainties. The first source of possible uncertainties is the assumptions in the compTT model 
(when used for the soft component): the electrons have a Maxwellian energy distribution, the 
scattering geometry is either spherical or disk-like, and the seed photons have a characteristic 
energy and are all originated from the center of the scattering medium. These assumptions, 
however, may not differ from the actual accretion disk structure significantly. First, because the 
inner region should be radiation pressure dominated, the disk temperature varies very slowly with 
radius (T(r) oc r -3 / 8 ), as opposed to a more rapid temperature variation in the gas pressure 
dominated disk (T(r) a r -3 / 4 ) (Liang 1999). Second, nearly all gravitational energy is released in 
a very small region around the inner disk boundary, thus the point and central source assumption 
is also not violated significantly. 

The second source of possible uncertainties is in the assumption that the soft component may 
be modeled as a disk-blackbody spectrum, in order for deriving the inner disk radius. Clearly the 
whole soft component cannot be described with the disk-blackbody model, due to the saturated 
Compton scattering in the disk region. However, since the main purpose of the model fitting is 
determining the flux in the soft component, and the disk-blackbody model does provide adequate 
fit to the high energy end of the compTT spectrum (as shown in Fig. 10), we believe it is a 
good approximation to use the diskbb model for the soft component. Therefore we conclude that 
color correction procedure we proposed is a reasonably reliable approach for determining the 
actual color correction factor from observation data. In fact the estimated color correction factors 
agree quite well with the possible values calculated from numerical analysis of the accretion disk 
structure (Merloni et al., 1999). 

The detailed physical processes occurring in the disk and corona are not addressed in this 
paper. We refer readers to another publication by some of us, in which we have identified a 
three-layered structure of the accretion disk, very similar to the solar atmosphere; between the 


commonly known cold disk and the hot corona, there also exists a warm layer, which is actually 
responsible for regulating the observed soft component, i.e., the color hardening process (Zhang et 
al., 2000). 


5. Conclusions 

We have re-analyzed some of the EXTE/PCA data of the microquasar GRO J1655-40 
and confirmed qualitatively the previously reported significant variations of the apparent inner 
accretion disk radius, before the color and radiative transfer corrections. However, our final 
results are quite different from the earlier ones, because we applied a physical model for the hard 
component ( compTT ), instead of the commonly employed phenomenological power- law model. 
The inner disk radius inferred directly from the spectral fitting varies dramatically between 
observations, if no Comptonization radiative transfer and color corrections are applied. After 
the Comptonization correction and color correction, the inner disk radius becomes quite stable. 
We therefore conclude that the inner disk radius in the microquasar GRO J1655-40 remains 
remarkably stable in a span of several years whenever the source is observed in a bright state; the 
corrected inner disk radius is apparently independent of the flux and spectral state of the source. 
(Note that our conclusion is irrelevant to the quiescent state of the source.) 

The stability of the inner disk boundary in GRO J1655-40 further supports that the compact 
object in the center of the accretion disk is a black hole, because the radius of the inner disk is 
expected to vary with mass accretion rate if the compact object is instead a neutron star, due to 
the interaction between the disk and the neutron star magnetosphere and/or the radiation drag of 
the surface X-ray emission of the neutron star. Using the ASCA data we are also able to derive the 
color correction factor from the data; the color correction factor changes significantly for different 
state of the source. Applying color corrections to the diskbb component derived with the ASCA 
data, we infer that the angular momentum of the black hole in GRO J1655-40 is around 0.8-1.0 
GAf/c 2 , consistent with some previous reports (Zhang et al., 1996, Zhang, Cui & Chen 1996; Cui, 
Zhang & Chen 1997; Kaaret ? ). 

The success of the simple Comptonization radiative transfer correction thus provides strong 
support to the Comptonization origin of the hard photons and indicates that the dominant source 
of seed photons to the Comptonization process in the corona is the soft photons from the disk. 

It is interesting to note that the outer truncation radius of the corona is quite large (800 r 9 ), in 
order for the inner accretion disk radius in the very-high/variable state to be consistent with that 
determined in the soft state. 

RXTE/PCA alone cannot be used to determine the color factor, due to the poor response 
of PCA below 2.5 keV. The inconsistency between the normalization parameter for the soft 
component between the RXTE/PCA and the ASCA/GIS instruments suggests that the absolute 
value of the diskbb normalization parameter (therefore the inner disk radius) cannot be determined 


with the RXTE/PCA data alone, again due to the poor low energy response of the RXTE/PCA 
instrument. 
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Appendix A. Radiative Transfer Correction to the Soft Component 

Neglecting some line and edge features in the X-ray spectrum of a black hole X-ray binary, 
the conventionally used two-component spectral model, i.e., multi-color disk blackbody component 
plus a power-law component with a high energy cutoff, is over-simplified in two aspects: 

First, the assumption that the power-law component extending straight to the low-energy 
limit of the spectrum is non-physical, since the power-law component is believed to be produced 
by inverse Compton scattering of low-energy photons by high energy electrons. The most natural 
source of the seed photons is the blackbody-like component produced from the accretion disk. 
Therefore a sharp low-energy cutoff or turn-over must be present in the power-law component 
in the low-energy region, since the blackbody-like soft component has a peak temperature, i.e., 
a turn-over in the low-energy end. Usually the observed peak temperature is between 2 to 5 
keV, thus requiring a sharp cut-off or turn-over of the power-law component roughly below 2 to 
5 keV. Neglecting this low-energy cutoff by applying a simple power-law spectrum model would 
under-estimate the flux in the soft blackbody-like component; the amount of under-estimation is 
apparently correlated with the shape and flux of the power-law component. 

Secondly, it has been a rather standard approach by many (including the authors of this 
paper) in the field to take the blackbody-like flux derived from the spectral model fitting as the 
true flux of the accretion disk. This ignores the radiative transfer process in the production of 
the power-law component. Since the power-law component is likely produced by scattering the 
soft photons in the original blackbody-like component of the accretion disk, each photon in the 
power-law component comes at the expense of a lost photon in the blackbody-like component, 
even if no absorption occurs in the Comptonization process. Therefore one must correct for this 
radiative process in order to infer the original flux of the blackbody-like component. Since the 
power-law production region cannot cover the entire accretion disk uniformly, the amount of 
correction to the blackbody-like component is thus energy dependent. Therefore ignoring this 
correction would result in under-estimation of the flux of the blackbody-like component. 

The above two problems tend to under-estimate the flux in the blackbody-like component, 
and are related to the co-existence of the hard X-ray power-law component. We thus expect that 
in the absence of any significant hard X-ray component (the conventionally called soft or high 
state), none of the above problems will cause any significant under-estimation of the blackbody-like 
flux. Indeed for all sources we know so far the inferred inner accretion disk radius is remarkably 
constant when no significant hard X-ray flux is observed (Sobczak et al., 1999a; 1999b). Therefore 
the inner accretion disk radius inferred when a system is in a high or soft state may indeed be the 
last stable orbit of the blackhole, although cares must be taken when the absolute value of the 
inner disk radius is estimated, since the absolute flux calibration of the instruments used must 
be known very well, as well as several correction factors must be applied (Zhang, Cui and Chen 
1997), even after the true luminosity of the accretion disk is derived from the observation. 

Ideally one would like to have a full physical model, in which the blackbody-like component 
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*phtons/cm 2 /s, 2.5-20 keV 
b 10“’® ergs/cm 2 /s, 2.5 - 20 keV 
c 2.5 - 20 keV 

d Electron number density in the inner disk radius. 

'Radiative transfer corrested normalization of dt»kbb component. Refer to Appendix A for detail. The outer radius of the corona is set to 800r e . 
The corresponding electon densities are listed in the last column of this table. 


Averaged ASM unit count rate for the day 



Fig. 1. — ASM light curve for the 1996-1997 outburst. The little vertical bars indicate the 
observations whose results are presented in this paper (Refer to Table 5). The dashed lines separate 
different states. Typical spectra for different state will be given for (a), (b), (c) and (d) (Refer to 
Fig. 2). 



Photons/cm 2 s keV ^ Photons/cm s keV 


±~t 


* 


960725 



10 

channel energy (keV) 
a) Very-high/Variable state (960725) 
970708 



10 

channel energy (keV) 

(c) High/Soft state (970708) 


> 960801 

QJ 

M 

w 10° 

w 

S io-' 

O 

^ 1 0 -2 
w 1 u 

5 io -3 

o 
A 

O- 2 
x 0 
-2 

10 


channel energy (keV) 
(b) Very-high/Variable state (960801) 
970825 



10 

channel energy (keV) 

(d) Low/Hard state (970825) 


Fig. 2. — Typical spectra in the three states. The upper panels are unfolded spectra, and the lower 
panels are residuals. The data from the four observations are fitted with the typical model we used: 
wabs(diskbb+compTT). The upper two spectra, which are all in Very-high/Variable state, are only 
several days apart, but they are fundamentally different; in (a) the soft component dominates, 
while in (b) the hard component dominates. The effect of radiative correction is very significant 
in the cases in which the hard component dominates, while less significant or neglegible for soft 
states. (Refer to Table 5 and Fig. 4 for the effect of radiative correction.) 
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Table 2: Fitting parameters from ASCA and XTE observation of 970226. (See Fig. 5) 



nH (10 22 cm~ 2 ) 

Tin (keV) 

Norm 

T 0 (keV) 

kT (keV) 

T 

ASCA alone 3 

0.63 

1.171 

1196 

- 

- 

- 

XTE alone b 

0.66 (fixed) 

1.125 

1817 

1.5 

49 

1.6 

Joint ASCA b 

0.66 

1.146 

1310 

0.14 

66 

1.6 

Joint XTE b 

0.66 

1.146 

1611 

Same as above 



°Hard component is not significant under 10 keV. 

6 A smeaxed Fe absorption edge component is included, with the energy fixed at 8 keV, width fixed at 7 keV. (Sobczak 
et al. } (1999)) 

and the power-law component are coupled together self-consistently, i.e., some of the soft photons 
produced by the accretion disk are up-scattered by the high energy electrons into the power-law 
component; the amount of the soft photons to be scattered (also in a frequency-dependent way) 
and the shape and flux of the power-law component depend upon the geometry and energy 
distribution of the corona. However, such a model is not currently available in the field, and may 
not be practical to make, since many degrees of freedom necessarily exist for this model, in the 
absence of the detailed properties of the hard X-ray producing corona. Here we propose for a 
semi-physical model for addressing the existing problems of the simple two-component model. 

Since the under-estimation of the flux of the blackbody-like component is most serious when 
the source is in a power-law bright state, whose power-law component usually exhibits a clear high 
energy cutoff above 100-300 keV, thermal Comptonization is the most likely mechanism for the 
power-law component production, although some non-thermal contribution, or a multi-temperature 
corona may be required in some cases. Because we are interested mainly in re-constructing the 
blackbody-like component reliably and the one-temperature thermal Comptonization model gives 
adequate description to the power-law component at low energies, we will only consider the 
one- temperature thermal Comptonization model for the power-law production in this paper. 

Titarchuk and Lyubaxskij (1995) have derived the analytical expression of the power-law 
component for the thermal Comptonization model in a self-consistent way, by assuming that 
the corona is in either a slab or spherical configuration with a uniform density and the soft 
photon source is located in the system of the corona. As will be discussed later, neither of 
these assumptions will likely introduce significant errors to the estimation of the blackbody-like 
component, as the final radiative transfer correction will at least partially compensate for the 
errors introduced by these assumptions. 

We therefore replace the power-law component in the simple two-component model by 
the Titarchuk and Lyubarskij (1995) model (the compTT model in the XSPEC package). The 
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Table 3. ASCA fitting parameters, derived color correction factors and the black hole spin. 


Date 


diskbb 


compTT 


/ b 

! 2 y/K^> 

h* 


nH a 

r in (keV) 

K bb nH a 

T 0 (keV) 

kT (keV) 

r 




94/08/ 23 c 



0.26 

0.20±8;81 

1.56 

8.1 


... 


Free nH 

94/09/27 

0.27 

1 4Q+ 0 * 005 
1 -^ y -0.004 

192±| 0.39 

o.2il°;°i 

1.21 

12.3 

2.62l°;l 3 3 

951^ 

0.88i°;° 9 9 

95/08/15 

0.69 

1 49+0-030 
i * 4O -0.020 

d 979t®g 0.53 

0.391°;°', 

1.20 

12.2 

i.36l°;°; 

58i; 

0.981°$ 

97/02/25 

0.63 

1 17+0 004 
A * A ' —0.005 

1196^21 0.54 

0.30l°;°2 

0.99 

13.4 


7211° 

0.94l°;| 3 3 

Fix nH e 

94/09/27 

0.27 

1 40+0.005 
J -’^ y -0.004 

192^2 0.27 

0-2641°°°; 

1.19 

13.0 

2.091°;°° 

60l 2 

0 Q7+0-03 
' -0.03 

95/08/15 

0.69 

1 49+0.030 
1 *^°-0.020 

0.69 

n oqo+0.005 

U.^O<£_ 0 .006 

1.17 

12.7 

1 OQ+0.00 

I.OO_Q 04 

nolj 

o.8oi°;°; 

97/02/25 

0.63 

1 17+0 004 
-0.005 

11 961 2 ? 0.63 

n ocn+0.008 

0.260_o.oo8 

0.99 

13.2 

1 67+ 0 * 05 
-0.05 

961° 

o.87l°;°5 


a in unit of 10 22 /cm 2 

b Color correction factor, / = T m /T 0 /2.7 

c Data from this observation cannot be fitted with diskbb model, so our method of deriving color correction 
factor cannot be applied to this data set. 

d Radiative transfer corrected. The original value is 786. 

e The nH is fixed to the value obtained from the fitting using wabs(diskbb) model. 
f Estimated black hole spin in unit of GM/c 2 . 


high-energy cutoff of the power-law component is also taken into account self-consistently in this 
model. The spectral model is thus 

F x = wab(diskbb + compTT) (2) 


Assuming that all photons in the power-law component are produced by Comptonization of 
the photons in the disk blackbody-like component by high energy electrons in a corona above the 
disk, it is in principle possible to derive the original emissivity of the accretion disk by carrying 
out a radiative transfer correction. It is already known that a simple uniform corona model cannot 
be used to fit the observed data; a density distribution in the corona is required, although the 
detailed geometry is still not clear yet (Hua, Kanzanas and Cui 1998). Here we assume that the 
corona has a spherical geometry with radius Ro and an electron density profile of: 

n e = n e fi{r g /R) p (r/i < R < Ro) (3) 

where r g = is the gravitational radius of the black hole, r h is the radius of the horizon of the 
black hole, p is index of the density profile and is taken to be unity for simplicity in this paper. 

In practice Rq may be determined approximately by examining the phase-lag spectrum of 
the light curve of a source. n e> o may in turn be determined by the optical depth inferred in the 
compTT model, or by comparing the blackbody-like and the power- law fluxes derived in the 
spectral fitting, as will be illustrated later. This assumed corona geometry is not very different 
from that inferred by Hua, Kanzanas and Cui (1998) for Cyg X-l, except that we also take into 
account the black hole in the center in a very simple way. 


We now consider a beam of photons emitted from a radius r < Ro on the disk surface and 
along the direction with an angle 9 from the disk surface. The optical depth for this beam can be 
described as: . 

r(r,8) = f n e a e dl (4) 

Jo 


where a e is the Thomson cross section, l 0 = -rcostf + \JRl~ r 2 sin 2 0, which is the length of 
the path when the beam intersects the boundary of the sphere for 0 < 8 < tt - arcsin if 
the beam reaches the black hole horizon bef ore intersects the boundary of the corona, i.e., for 
7 r - arcsin ^ <8 < ir, then l 0 = -rcosd- yV); - r 2 sin 2 0. 


Substitute Eq. (3) and R = \/r 2 + 1 2 + 2rl cos 8 into Eq. (4), we obtain: 

(0 < 8 < 7T — arcsin ^ ) 
(tt — arcsin ^ < 8 < n) 


, n\ 1 Rc + y/Rf^-^e 

n(r, 8) = n efi a e r g In 

l n Wr»-rW 
1 m r(l4-cos 9) 


T2(r,0) = n e fla e r g ' 


( 5 ) 


The average probability P(r), which describes the probability that a disk photon is not 
scattered by an electron in the corona, can be calculated from Ti(r, 9) and T 2 (r, 6): 

^-arcsin ^ e _ riM ) sin ^ + ^ r „ e^M) sin 8d9 

Jq sin 8d9 


{Tin < T < Ro) (6) 


where r>in is the inner radius of the accretion disk. 

If the photon is emitted from the disk surface outside the corona sphere, the average 
non-scattering probability P{r) can be expressed as: 


fr* 1 ™ * sin Odd + e- Ti(T ' d) sin (6)d6 + ^ e~^ s in(0)d0 

P ^ = " Jo" sin 9d§ 

where r 3 (r, 0) and r 4 (r,0) is described as 


(r>R o) 
(7) 


T 3 (r,0) = n e ,oa e rg\n 


t 4 (r, 6) = n e fia e r g In 



(7T — arcsin ^ < 9 < n — arcsin 
(7r — arcsin ^ < 6 < 7r) 


The equivalent optical depth f(r) can then be calculated from 


( 8 ) 


f (r) = - In P{r) 


(9) 


where P(r) may be calculated from Eqs. (6) and (7) for r < Ro and r > Ro, respectively. In Fig. 7, 
the equivalent optical depth f (r) is depicted as a function of r / Ro for different combinations of 
the Tin and rh corresponding to different values of the black hole spin angular momentum (see 
Zhang, Cui and Chen 1997 for details). The slight turn over near the inner disk boundary is due 
to the existence of the black hole horizon, although no other general relativity effects are taken 
into account for the Compton scattering process near the black hole horizon. 


In the multi-color disk blackbody model, the local disk radiation is assumed to be blackbody, 
i.e., its surface brightness follows the Planck distribution B(T(r), u) = i}> an ^ 

T(r) = T in (r/r tn y 3 ^. Therefore the average non-scattering probability for a photon with a 
frequency v can be calculated as 


_ J” 27rrE(r,^)P(r)dr 
S£2*rB{r t v)dr 


( 10 ) 


The equivalent optical depth for a photon with a frequency v is thus 


f{u) = -\nP{v) 


( 11 ) 


If both Ro and n e , 0 are known, then P( v) and thus t(v) can be calculated. It is easy 
to show that r(r) ~ n e fiQ, e r g ln(7?o/ 7 'm)i if r ~ r» n and iZo i"m. It then follows that 
f{v) ~ n e fla e r g ln(iZo/rj n ), if hv ~ kT in . This also implies that n efi and Ro are coupled together 
in n e fia e r g ln(Ro/ri n ). In Fig. 8, the equivalent optical depth t{v) is depicted as a function of 
the emitted photon energy, for different combinations of n e> o and Ro, which are coupled together 


by n efi a e r g \n{Ro/r in ) = constant ; the initial value of the constant is determined by taking 
Rq = 10 5 , n e> o = 10 17 /cm 3 (all radius are in units of the gravitational radius of the black hole, i.e., 
r g — GM/c 2 ); kTi n = 1.0 keV. The left figure corresponding to the Schwarzschild geometry, i.e., 
ri n = 6 and = 2; the right one is for a near-extremal Kerr black hole with rj n = 2 and r* = 1. 

Once P(v ) is known, we can derive the initial disk spectrum Fi(v) from the observed disk 
spectrum F 0 (u) (the diskbb component in Eq. (2)): 

Fi{v) = (12) 

COS l 

where i is the inclination angle of the accretion disk. 

As we mentioned above, Rq may be estimated by examining the phase-lag spectrum of its 
light curve. On the other hand, the value of n C) o may be calculated from the following radiative 
transfer equation: 

[ F(compTT, v)dv = f F 0 (diskbb, i/)(e f 1 M - 1 )dv (13) 

J COS % J 

where F(compTT, u) and F 0 (diskbb,v) are the observed power-law and blackbody-like photon 
spectra, respectively. 

Finally we can use the corrected accretion disk emission spectrum to estimate the inner 
accretion disk radius. First we can obtain the total blackbody-like flux: 

F(BB) = J huF 0 (diskbb,u)e f ^du (14) 

Then we can apply the Eq. (3) in Zhang, Cui and Chen (1997) to estimate the inner accretion 
disk radius, i.e., the radius of the last stable orbit of the black hole. 

It should be noted that the escaped spectrum also has a slightly different shape than the 
original spectrum, due to the frequency dependency of the equivalent optical depth f(u). In Fig. 9, 
we show the escaped spectrum, i.e., the observed one, corresponding to different configurations 
of the corona, for a Schwarzschild black hole ( left figure: rj n = 6 and = 2) and a Kerr black 
hole ( right figure: rj„ = 2 and r^ = 1). The values of n ej o and Rq are not allowed to change 
independently, i.e., are coupled together by n e Q ln(fZo/rj n ) = constant (the initial value of the 
constant is determined by taking Rq = 10 5 , n e> o = 3 x 10 17 /cm 3 ); only values of Rq Eire marked 
in the figures. The original disk spectrum without any scattering is also shown for comparison 
(kTi n = 1.0 keV). It is clear that the difference of the spectral shape from the original spectrum is 
in general very small, especially at energies > fcTj n , unless when Rq ~ 10 2 . For the Kerr black hole, 
the escaped spectral shape is closer to the original one for the same corona configuration; however, 
the same spectral shape may be produced for both types of black holes, by taking different corona 
configurations. 


Appendix B. The Color Correction 
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Even after the soft component is properly extracted from the observation data, proper color 
correction to the component is still required, since the emerged spectrum from the disk is actually 
not purely blackbody, due to electron scattering in the disk. In the optically thick disk the electron 
scattering is saturated (optical depth significantly greater than unity), thus the escaped spectrum 
may still be described as a multi-color blackbody model, in which the local spectrum is described 
by a diluted blackbody spectrum: 

B'(hu, T'(r)) = fcoiT{r)) (15) 

■'col 

where B is the Planck function and fcol is the color correction factor Therefore the inferred disk 
radius is proportional to Previously the value of fcol has been estimated using theoretical 
calculations and assumed a constant in most cases. However, since the real value of the color 
correction depends on the property of the electron gas in the disk, in principle the color correction 
factor may not be a constant for a source in different spectral and luminosity states. 

We may also model the soft component with the compTT model, from which we can obtain 
directly the temperature (To) of the seed photons. If the seed photon spectrum is diskbb, instead 
of blackbody assumed in the compTT model, the real inner disk temperature of the disk should 
be kTo * 2.7. Thus according to the definition of the color correction, f^i = T* n /2. 7/cTo, in which 
Ti n is obtained by fitting the same soft component with the diskbb model. The consistency of this 
approach may be examined by comparing Tj n and T e , which is the electron temperature inferred 
in the compTT model fitting; we should have Tj n ~ T e if the observed soft component is indeed 
produced by saturated thermal Comptonization in the disk. In Fig. 10 we show one simulation for 
the ASCA GIS instrument, in which the input spectrum is generated from the compTT model 
(parameters shown in the figure), and the data are then fitted with both compTT and diskbb 
models. It is clear that the diskbb model also provides adequate fit to the data, and Tj n « T e 

This offers a simple way for estimating the color factor from the observation data, thus 
avoiding many uncertainties in theoretical calculations and the assumption that the color correction 
factor is universally the same for all sources in all observations. There is, however, one restriction 
for applying this color correction method: the low energy response of the instrument must be 
good enough to allow kTo to be determined accurately. We have carried a series of simulations 
for evaluating capability of the existing instruments. We found that the ASCA instruments are 
the best suitable for this purpose, though the SAX instruments are comparable in performance. 
However, the RXTE instruments do not have the adequate low energy response for estimating the 
color correction factor this way. 

We thus analyzed some ASCA archival data on GRO J1655-40 when the source was active. 
The results are summarized in figures 5-8. It can be seen that the corrected inner disk radius for 
all observations is consistent to a constant, despite the large spectral and flux variations between 
these observations over several years. This demonstrates that this technique in estimating the 
color correction factor is quite reasonable and that the inner accretion radius is indeed very 


stable. We argue this as the most convincing evidence for the stability of the inner accretion disk 
boundary in a black hole X-ray binary. 
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Fig. 3. — Comparison between power-law and compTT model. Notice the cut-off at low energy end 
of the compTT model and the slightly underestimation of the diskbb normalization. 







MJD - 50200 

Fig. 4. — The main results of our work. It can be seen that the observed soft component 
normalization varies dramatically in the very- high/ variable state. Also there is a significant anti- 
correlation between the hard photon flux and the soft component normalization, which indicates 
that they must be of the same origin and provides evidence for the necessity of Comptonization 
radiative correction. The bottom panel shows the soft component normalization after the correction, 
which is much more stable. Photon flux in this figure is in 2.5 - 20 keV. 
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Fig. 5. — Joint fitting of ASCA and XTE data from the observation taken in 970226. The 
ASCA/GIS2 data in 0.7 - 10 keV and XTE/PCA data in 2.5 - 20 keV are used. There exists 
a significant difference between the ASCA and XTE normalizations (Refer to Table 2). 
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Fig. 6. — Three ASCA observations taken in 970226, 950815 and 940927 (from top to bottom), 
when the source is in soft state and the hard component is negligible, therefore the Compton 
radiative correction is not necessary. It can be seen that the spectral shape and normalization 
varies significantly. However, the derived color correction factor is not constant, the inferred since 
the but the inner disk radius (after color correction) remains rather stable (Refer to Table 5). 
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Fig. 7. — The equivalent optical depth f(r) is depicted as a function of r/Ro for different 
combinations of the n„ and r/i corresponding to different values of the black hole spin angular 
momentum (see Zhang, Cui and Chen 1997 for details). The slight turn over near the inner disk 
boundary is due to the existence of the black hole horizon, although no other general relativity 
effects are taken into account for the Compton scattering process near the black hole horizon. 
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Fig. 8. — The equivalent optical depth t (v) depicted as a function of the emitted photon energy, for 
different combinations of n C) o and Rq, which are coupled together by n e fla e r g ln(/?o/r m ) = constant ; 
the initial value of the constant is determined by taking Rq — 10 5 , n e ,o — 10 17 /cm 3 (all radius are 
in units of the gravitational radius of the black hole, i.e., r g = GMjc 2 ); fcT< n = 1.0 keV. The left 
figure corresponding to the Schwarzschild geometry, i.e., rj„ = 6 and r/, = 2; the right one is for a 
near-extremal Kerr black hole with = 2 and = 1. 
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Pig 9 — The escaped spectrum, i.e., corresponding to different configurations of the corona, for 
a Schwarzschild black hole ( left figure : r in = 6 and r h = 2) and a Kerr black hole ( right figure: 
n n = 2 and r h = 1). The values of n e ,o and Ro are not allowed to change independently, i.e., are 
coupled together by n efi ln(Ro/r in ) = constant (the initial value of the constant is determined by 
taking Ro = 10 5 , n e>0 = 3 x 10 17 /cm 3 ); only values of Rq are marked in the figures. The original 
disk spectrum without any scattering is also shown for comparison ( kT in = 1.0 keV). It is clear that 
the difference of the spectral shape from the original spectrum is in general very small, especially 
at energies > kT in , unless when Rq ~ 10 2 . 
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Fig. 10. — Simulation for fitting thermal Compton spectrum with diskbb model in XSPEC. The 
simulated data are produced with the compTT model. The data are then fit with both compTT 
(left) and diskbb (right) models, both producing adequate fitting. Data parameters: T 0 = 0.30, T e = 
l.OkeV, r = 10; Fitting parameters with compTT model: To = 0.29, T e = 0.94keV,r = 11.0; Fitting 
parameters with diskbb model: T\ n = 0.96keV. 






